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FORCES IN NONLINEAR MEDIA 


Mordehai Milgrom 

Department of Condensed Matter Physics, Weizmann Institute, Rehovot Israel 

I investigate the properties of forces on bodies in theories governed by the generalized Poisson 
equation V • [/r(|V</9|/ao)V(^] oc Gp, for the potential ip produced by a distribution of sources p. 
This equation describes, inter alia, media with a response coefficient, p, that depends on the field 
strength, such as in nonlinear, dielectric, or diamagnetic, media; nonlinear transport problems with 
field-strength dependent conductivity or diffusion coefficient; nonlinear electrostatics, as in the Born- 
Infeld theory; certain stationary potential flows in compressible fluids, in which case the forces act 
on sources or obstacles in the flow. The expressions for the force on a point charge is derived exactly 
for the limits of very low and very high charge. The force on an arbitrary body in an external field 
of asymptotically constant gradient, —go, is shown to be F = Qgo, where Q is the total effective 
charge of the body. The corollary Q = 0=>F = 0isa generalization of d’Aembert’s paradox. I 
show that for G > 0 (as in Newtonian gravity) two point charges of the same (opposite) sign still 
attract (repel). The opposite is true for G < 0. I discuss the generalization of this to extended 
bodies, and derive virial relations. 

PACS numbers: 


I. INTRODUCTION 

The Poisson equation, which governs so many physical processes, has the nonlinear generalization 

V -[pilVip\/ao)Vif] =anGp, ( 1 . 1 ) 

by which the source distribution p(r), in U-dimensional Euclidean space, gives rise to a potential held tp. Here, oq is 
a constant with the dimensions of Vcp, an = 2(7r)^^^/r(Il/2) is the H-dimensional complete solid angle, introduced 
here for convenience, and G is a coupling constant. As I will show, for G > 0, a point, test charge is attracted to a 
(hnite) point charge of the same sign (as in gravity), while for G < 0 it is repelled. 

Equation(^r^) describes a variety of physical problems; some examples are: 

(i) Nonlinear dielectric, and diamagnetic, media; p is then the dielectric or diamagnetic coefficient, which depends 
on the held strength (here G < 0). 

(ii) Problems of nonlinear electric-current hows in systems with held-dependent conductivity (nonlinear current- 
voltage relation), and nonlinear dihusion problems; p{w) is the transport coefficient. 

(iii) Stationary, subsonic, potential-how problems of non-viscous huids with a barotropic equation of state p = p{g) 

(j> is the pressure, g the density). The stationary Euler equation is integrated into Bernoulli’s equation f{g) = 
— + const., where f'{g) = g~^p'{g) = <?{g)/ P: with c the speed of sound; / is thus increasing with g, and p is a 

function of |u| = |V:/?|. The stationary continuity equation then gives V-[e(|V:p|)V:/3] = s(r), with s the source density 
(see e.g. |Q for the ideal-gas case). This is eq.(pTl|) with p = g, and G = > 0. For example, in a huid with an 

equation of state of the form p = ag'^ (a > 0, 7 > 1), we have g{u) = p(0)[l — {u/uo)^Y^G-D ^ ^ith itg = 2co/(7 — 1), 
and Co is the speed of sound at w = 0. Subsonicity requires (u/uq)^ < (t ~ l)/(7 + !)■ 

If we, formally, consider a stationary how problem in a medium with negative compressibility, < 0, ellipticity is 
maintained for any value of 'S/ip. For example, for a medium with a constant < 0, g(u) = g{0)exp{u'^/2\c\'^). 

(iv) Nonlinear (vacuum) electrostatics as formulated e.g. in the Born-Infeld nonlinear electromagnetism, which also 

appears in ehective Lagrangians resulting from string theory (see review and references in |Q] 0). In the original, 
electrostatic Born-Infeld theory p{w) oc (1 — and G < 0. 

(v) A formulation of an alternative nonrelativistic gravity to replace the dark-matter hypothesis in galactic systems 
1^. Here p{w) « w for w <C 1, and /i ~ 1 for w » 1 (G > 0). 

(vi) Equation( |l.l[) was used in as an effective-action approximation to Abelianized QCD. 

(vii) Area (volume) minimization problems, such as the determination of the shape of a soap him with a dictated 

boundary (see e.g. [^): If Xn+i = p{xi,Xn) describes a D-dimensional hypersurface embedded in D-l-1 dimensional 
Euclidean spa ce with Cartesian coordinates Xi ,..., Xo+i, the volume element on the surface is dv = + d^r. 

Then, eq. 0 ) describes the problem of the minimization of the volume of the surface. The sources may describe a 
force density on the hypersurface acting in the direction Xn+i- In this problem p{w) oc (1 -I- and G > 0. 

Born-Infeld electrostatics is the same as the area-extremization problem for a time surface embedded in Minkowski 
space-time. 
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Much has been said in the mathematical literature on properties of the potentials that solve eq. ( 0 ) (see e.g. [^). 
But, to my knowledge, very little has been said about forces on bodies in such theories. The forces can be written as 
certain integrals of tp, and are of obvious relevance in the physics context. 

Many of the familiar and intuitive properties of the linear theory are lost in the nonlinear case because the potential, 
and forces, is not the sum of the contributions of the sub-systems. For example, the force on a point charge is no 
more proportional to the charge, does not reverse direction when the charge of the body reverses sign, etc.. And, 
Earnshaw’s theorem [|j no longer holds. Under some circumstances it is possible to suspend stably static charged 
bodies in a static field. This last aspect is treated in detail in 0. 

The choice fi{w) = w^~‘^ is special in that the theory is then conformally invariant, and lands itself to many 
analytical developments. This case is describe d in detail in Q. 

After discussing some general aspects of eq. (0 in section II, I take up the main subject concerning the properties 
of forces on bodies in the nonlinear theories: general properties in sections III, and forces on point charges in section 
IV. I conclude in section V with some examples of applications. 


II. GENERAL PROPERTIES 


A. Preliminaries 


The field equation ( 0 ) is derivable from the action functional 

S = S, + Sf = -J pip d^r - J .F[(Vy))7ag] d^r. (2.1) 

The function p{w) in eq. (0 is given by 

’ y = w‘^. ( 2 . 2 ) 

ay 

For equilibrium problems, such as (i), (iv), (v), and (vii) above, the quantity E = —S may be identified with the 
energy of a charge configuration p when S is finite. I shall only be interested in differences in the energy between 
configurations with the same total charge, and this is finite, in general, even if the expression for E diverges. In the 
case of non-equilibrium, stationary transport problems, such as nonlinear diffusion, heat transfer, etc., E is not an 
energy but is related to the entropy-generation rate. 

For convenience, the free, additive constant in T is chosen such that 1F(0) = 0. This makes the field-action density 
vanish at r —> cx) when V:/? —> 0 there, and improves convergence. I exclude from discussion theories for which T 
diverges at zero argument. I also assume p{w) > 0 except, possibly, at ru = 0 where p may vanish. Thus T(yj) is an 
increasing function, and is positive for y > 0. (I assume that E'{y) does not vanish for y > 0; it then has a uniform 
sign which can always be taken as positive by adjusting the sign of G.) 

In the area-minimization case Eiy) = (1 -I- y)^/^ — 1. In the flow problem Tiy) is essentially the equation of state 
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since the pressure is given by p = p(u = 0) — ^T[u'^[q)/u^, with the additive constant chosen so that .F(O) = 0 

[p < p(o)]. 

All that I say below is straightforwardly generalized to the case where IF, and thus p, depends explicitly on r. 
For instance, in the flow problem this happens when the fluid is coupled to some potential field "0 via g{r)ip{r) for 
which the Bernoulli equation becomes f{g) = — '*A(r) -I- const.. To avoid encumbrance I assume all along, unless 

otherwise said, that there is no such explicit r dependence. 

When the charges are not held fixed but move under the influence of the ip field, and their dynamics is of interest, 
we add to the action the kinetic term for the charges 


S^ = 


d^'rg.mV 


(2.3) 


where Pm is the mass density of the particles contributing p to the charge density. Extremizing S Sp, with respect 
to particle coordinates give the usual Euler equation g^^ = —pVtp. 

The solution of eq.(l.l) is unique inside a volume V when one dictates on its boundary the value of ip or that of 
p{\Vip\/aQ)dnip (or a combination thereof) (i9„(p is the normal component of ip) provided the function i>{w) = wp{w) 
is an increasing function (see e.g. |^] for a proof). This monotonicity of ^(i(;)-which I shall assume all along-is 
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tantamount to /t = din ^{w)dln w > —1. This is also the condition for the ellipticity of eg. ( ^.11 ). In the case of a 
stationary flow we have fi = —u^/c^; the ellipticity condition is then equivalent to subsonicity of the flow. 

I now show that from the ellipticity condition and the choice 1F(0) = 0 follows that the logarithmic derivative of JF, 
T{y) = yT'{y)/T[y), satisfies !F[y) >1/2 for all y > 0 for which ellipticity obtains. I use this inequality repeatedly 
in what follows. Define 


x{w) = — 1] = 2y{w)w'^ — (2.4) 

(where y = w^). First note that x'(rc) = 2wy{l + fi), so x' > 0 for w > 0. Since x(0) = Oj x(w^) > 0 for ic > 0. 
Thus T{y) > 1/2 for all y > 0 [as T{y) > 0 for y > 0]. (The condition .F(O) = 0, on which the derivation of JF > 1/2 
depends, does indeed enter in the cases where this inequality is used below.) 

It is useful to write our theory for a general curved space whose metric is y^j {g'^ its inverse, and g = \det{gij)\). 
The covariant form of the action is then 


S = - 


J d^r - J y^/^J^[((^, V,,)/ao] d^r, 


(2.5) 


with if/ = g'^ifj = difldx"^), and p* = g is a scalar under general coordinate transformations. Repeated 
indices are summed over. The covariant form of the held equation is 


[F(F,V,i/ao)F,1;». = aoGp*, 


( 2 . 6 ) 


where a semi-colon signifies a covariant derivative. The covariant divergence appearing in eq.(^) is given in terms 
of the normal divergence of a vector as = g~^^^ig^^‘^y^),k- So, using usual derivatives instead we have 


= ttoGp. 


(2.7) 


From the covariant action we can derive the the field stress tensor (the energy-momentum tensor when working 
in Lorentzian space-time). This is the functional derivative of the held action with respect to the metric: under a 
variation Sgtj 


<5^/ = ^ / d^r. ( 2 . 8 ) 

In the Euclidean case, on which I concentrate hereafter, we And 

2 2 

V = ® W) = ® (2-9) 

2iO!.£)(jr 

where e = Vy>/|V(^| is a unit vector along V 93 . The trace of V is —{a^/2ariG)!F{D — 2.F). The held direction e is 
an eigenvector of V with eigenvalue —(aQ/2aDG).F(l — 2jF). The inequality T > 1/2 implies that this eigenvalue is 
always positive: there is always tension along the held lines. All other eigenvalues are equal, and negative. 

For solutions of the held equation the divergence of "Pij, which measures the rate of change of the momentum 
density, is given by 

V-P = pV(p. (2.10) 

This conservation law can be derived directly from the Euclidean action dH) and follows from its translation invari¬ 
ance: under infinitesimal translations r ^ r -|- a, yj(r) ^ (p{r) + (a • V)(^, etc. S^S + a.-{V-V — pVp). 

The field equation can also be written as 


= aoGp, ( 2 . 11 ) 

where 

A = (1-I-ye (g) e), (2-12) 

with /t the logarithmic derivative of y, and e (g e is the metrix whose {i,j) element is CiCj (all dependent on Vy). 
Since y > —1 is positive definite. If y(0) = 0, points where Vp = 0 need special treatment which I do not go into 
here. 
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If we make a small change 5p in p, the field equation can be linearized in the small change C in tp to read 

V ■[pA-^(\=ar>G5p, (2.13) 

This is the same as the equation for the electrostatic potential produced by the density 5p in a linear dielectric medium 
with a position-dependent, anisotropic dielectric constant pA. 

A variation of p in the volume V gives rise to the a variation in S 

6S =[ Sp{V ■ [p{\Vip\/ao)V(p] - aoGp} d^r - ^ [ pSpVip-ds, (2-14) 

oidG Jy aoG Jj, 

where S is the boundary of V. For potentials that solve the field equation, and hence nullify the first term in eq.( ^.14 ), 
we have 


6S = - — f p SipVip ■ ds. (2-15) 

J s. 

We can obtain useful integral constraints on such solutions -such as conservation laws and virial relations-by consid¬ 
ering specific variations that do not nullify the surface term. Some examples are given in Appendix A. 

The second-order change in the action (energy) is 

5‘^E = -6^S = [ p^5p ■ A ■ VSp d^r. (2.16) 

2aD<b Jy 


The ellipticity condition makes the integral positive (when Vtp ^ 0), and thus a solution of the field equation is a 
minimum of the energy for G > 0, and a maximum for G < 0. 

Using the integral relation (A.l) derived in Appendix A, which holds for solutions of the field equation, to eliminate 


the explicit dependence of E on the sources we get 


E = -S = - 


2aoG . 


T[{Vpf/al]{2T-l) d-r. 


(2.17) 


So, in light of the above inequality for E, E is positive for G < 0 (as in electrostatics). 

The field equation enjoys a certain scaling property Q in that if p{v) and p{v) are a consistent pair then so is 

(pA(r) = Av3(A-'r), pxiy) = A'V(A^'r), (2.18) 

with the appropriately scaled boundary conditions. Charges then scale as q\ = 


B. Asymptotic behavior of the potential 

When the medium can be considered infinite, a common choice of boundary condition, describing an isolated 
system, is > 0 at infinity. In fact, if the potential can be assumed to become spherical at infinity, this boundary 
requirement sometimes follows from potential equation itself through Gauss’s theorem. The behavior of p{w) near 
w = 0 is then relevant. For concreteness, I assume in the rest of the paper that T, and thus /i, approaches a power 
of the argument near 0 . 


p{w)^W^-, (/ 3 „>_ 1 ). (2.19) 

This choice does not cover the case in |Q], where p is logarithmic at small values of the argument. 

If the sources p are contained within a finite volume, and the total charge, Q, does not vanish, the field becomes 
radial at infinity, and, applying Gauss’s theorem to the field equation for a sphere of a large radius r, we find 
asymptotically 


V(/?« s(QG)|G|'>''>|QP'>r-^°<°-'>n. 

Here G = Goq”, 70 = 1/(1 -|- /3o), s(a) = sign{a), and n = r/|r|. 


( 2 . 20 ) 
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The value /3o = -D — 2 is a limiting one. For higher values of /3o, ip diverges like a power of r at large r (when Q ^ 0) 
with Vip still vanishing there. For the limiting case, ip is logarithmic at infinity, and V(^ oc rr“^. In this case the 
action is infinite when calculated for the whole space. The revised theory of Newtonian gravity discussed in [0 is of 
this limiting-power type, and so are the class of conformally invariant theories with fi = D — 2 discussed in ^ . In 
what follows I assume Po < D — 2 unless otherwise stated (as in the discussion of one-dimensional systems). 

For Q = 0 I do not have a general expression for the asymptotic behavior of the held. Because of the nonlinearity, 
multipoles are no more very relevant. For (3o = D — 2 the asymptotic form can be obtained as follows. In the 
exact-power-law theory with (3 = D — 2 one can use the conformal invariance. Take the origin at a point outside 
charges, and make a conformal transformation r — > a^r/r^, where a is the radius of the rehection sphere. If ip{r) is 
the solution of the original problem then (p{r) = ip[a3x jr^') is the solution of the new problem with the transformed 
charge distribution /5(r) = which also has a vanishing total charge (see for more details). So the 

asymptotic behavior of ip is obtained from the behavior of p at the origin (where the charge density is 0). I assume 
that (p is analytic there, so its dominant behavior is, generically, p w o“^K • r, where K is some constant vector. This 
gives for the generic, asymptotic behavior of p 


(p Ri K • r/r^. 


( 2 . 21 ) 


K might be viewed as the asymptotic-behavior dipole, but it is not proportional to the dipole of the charge distribution. 
If K = 0 the asymptotic behavior is, more generally, of the form 

( 2 . 22 ) 


where a totally symmetric constant tensor, is the first non-vanishing Taylor coefficient in the expansion 

of p at the origin. (I am discussing only the leading behavior, not expansion terms. Because of the nonlinearity the 
higher-order behavior depends on the leading one.) Basically, I use the conformal invariance to argue that p has to 
be analytic in r/r^ at infinity, and hence to greatly constrain its form there. The Ks have further to satisfy algebraic 
relations insuring that expression(2.22) satishes the vacuum held equation. For n=l there are no extra relations on 
K. For n = 2 the algebraic relation in the linear {D = 2) case is the usual tracelessness requirement Ku = 0. For 
D > 2 the condition on the matrix 1C, the elements of which are Kij, is {rlCr)^^~‘^'>/‘^r[{D — 2)K? -I- Trace{JC)K?]Y = 0 
for all vectors r. This, for the symmetric /C, can be shown to imply /C = 0, so there is no dominant n = 2 behavior 
for > 2. Now consider a general theory with (3o = D — 2. Let p be the solution for a conhned charge distribution 
with vanishing total charge. Dehne 


p* = {aoG)-^y ■[p{\yp\lao)^p\, 


(2.23) 


where p, is the exact D — 2 power. Because asymptotically (where ^ 0) we have p —>■ p, we get from Gauss’s 
theorem that p* also has a vanishing total charge. And, if p approaches its power-low behavior fast enough, p* will 
be well bounded. But from eq.( 2.23| ) v? is a solution for p* i n the theory with the exact power-law p and so, from the 
arguments above, must also have the asymptotic form( p)^ )(2 .22| ). 

Outside a spherical distribution of zero total charge the held vanishes for all theories. 

Also of interest is the boundary condition Vp ^ —go for r ^ oo (go is a constant vector). It pertains e.g. to a 
system of charges in an external electric held, to magnetized or superconducting bodies in an ext ernal magnetic held, 
or to obstacles and sources in an asymptotically uniform how. For the asymptotic held use eq.(2.13) to linearize in 
VC = Vtp-bgo: 


V-[VC + Aeo(eo-VC)] =0, (2.24) 

and solve with VC —> 0 at inhnity (bq = go/|go|)- Taking eo to be the a;i-direction we see that, asymptotically, C 
satishes the Laplace equation in the coordinates r' = [(1 -|- p)~^Gxi,x 2 , ■■■, Xd])- Thus, C has the standard multipolar 
asymptotic expansion in r'. When the total charge, Q, is hnite the dominant behavior is C oc {D > 2). 


III. FORCES ON BODIES 

A body is an isolated region, v, where the held is externally disturbed in one way or another. For example, a body 
may be dehned by some rigid distribution of charges in u, or by dictating the potential or its gradient on the boundary 
of V. The body may also be dehned as a “p inclusion”: dictating in u a p{w) that is different form the ambient one. 
Examples of bodies dehned by boundary conditions are: a rigid body in the how problem and a superconducting 
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body in a magnetic field, for both of which the normal component of vanishes at the surface; a conducting 
(equipotential) surface in an electric field; and a restricting boundary in the volume-minimization problem, for which 
the potential is dictated. Examples of bodies defined by fi inclusions are a dielectric inclusion, or, in the flow problem, 
a region of space where a fluid is subject to an external potential ^(r) that couples to the fluid density. There is then 
an extra term, ^(r), in the Bernoulli equation, and g{u^) is then also a function of r in u through ip. 

Consider then a background having some ambient, position-independent /i. Within this background are embedded 
disjoint bodies of different types as defined above. Once the solution, (^, it gotten for the system one can replace the 
bodies by effective charge distributions that give the same potential field outside the bodies. The effective charge 
density is simply 


p*{r) = (a,G)-iV • [/x(|V(^|/ao)V¥>], (3.1) 

with the ambient /i used everywhere. If a body is defined by boundary conditions we continue the solution inside the 
body by solving the field equation with the ambient p, and imposing the same boundary conditions for the internal 
solution. Obviously, the effective charges given in th is w ay appears as an extra charge only on the bodies because, by 
the field equation, the right-hand side of expression ( |3.l| ) gives the true charges outside bodies. For example, a body 
that is defined by dictating p on its surface is replaced by a p* that constitutes a surface charge; a body defined by 
dictating pdn^p on the surface is replaced by a surface dipole layer (whose total charge obviously vanishes). A body 
defined by a /i inclusion is replaced by an effective charge distribution whos e to tal effective charge is the same as the 
actual charge on the body. This is because Gauss’s theorem, applied to eq. gives the total effective charge as a 
surface integral over pVp on a surface surrounding the body, but, from the field equation, this also equals the actual 
total charge of the body. So the total effective charge of a /i inclusion always vanishes if there are no actual sources 
within it. 

The force on a body can be defined in several equivalent ways. For example, it may be taken as the gradient of 
the total energy under translation of the body: First replace all bodies in the system by their effective charges. Then 
translate the body in question, rigidly and infinitesimally, by Sa, keeping all the other charges (including effective 
ones) fixed. For the energy increment 6E = —6a • F„, we identify F„ as the force on the body. Under the translation, 
the effective charge distribution, p, changes by Sp* = -Sa - Vp* in v, and Sp* = 0 outside. Since p is an extremum of 
S, the change in E (= —S) can be calculated as if only p* has changed {5p = 0 at infinity). Thus 


(5a • Fu = —5E = — pSp* d^r = 6a- Vp* p d°r = —6a ■ / p*V(p d°r, 


(3.2) 


where the last equality is obtained by integrating by parts. Thus, 


F„ = - / p*Vp d^r = / pVp* d^r. 


(3.3) 


The first equality may also serve as a definition of the force: it says that the force is the sum of forces on the elements 
of the body, each given by —p*Vp d’^r. But, unlike the linear case, here p cannot be taken as that due to the rest of 
the system (excluding the body from the sources). In the linear case, but not here, the potential can be written as 
the sum of the contribution of the body and of the rest of the system. Since a body does not exert a force on itself 
(this is also true in the nonlinear case, see below) the contribution of the self field drops from the expression for the 

force. _ ^ ^ „ „ 

A third way-perhaps the most useful-to write the force uses eq. (|2.10|) (p*V(p = V • P) and expression ( ^.9[) for V 
to get 


F„ = - / P • ds = 


2aoG 


[Eds- / 2EE\Vp\-‘^VpVp-ds]. 


(3.4) 


The integration is done over any closed surface, a, that surrounds the body and excludes all other sources and bodies 
(compare with the expression of the force as a surface integral in |^). This expression has the advantage that it does 
not require the effective charge distribution of the body, and employs only the field outside the body. 

The surface integral in eq.(3.4) vanishes automatically for the surface at infinity in theories with Pg < D — 2, so 
the total force on a whole system vanishes, as expected from translational invariance. In other theories (e.g. in one 
dimension-see below) the boundary conditions for an isolated body must insure this. 

For nonlinear dielectrics, diamagnetics, etc. the above definition of the force coincides with the usual expression 
for the electromagnetic force. In the flow problem F^ is the mechanical force acting on a region containing sources (p 
replaced by s) because the integral in expression (|3.3|) (f su) is the rate at which the sources in the volume impart 
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momentum to the flow. For a rigid obstacle standing in the flow, expression (3.4) gives the mechanical force the flow 
exerts on the obstacle because the second integral vanishes as Vyi • ds vanishes on the surface (the flow is parallel to 
the surface) and we are left with Fy = — f pds. The same is true for a type-I superconducting body in a nonlinear 
magnetized medium: due to the Meissner effect Vtp • ds = 0 on the surface of the body. 

In the volume-minimization problem the F is the actual lateral force (in the xi,...,Xd plane) on the volume in 
question due to the uneven tension. 

In the case of non-equilibrium, stationary-transport problems, such as nonlinear diffusion, the action functional 
is not an energy, and F is not a force, but measures the gradient in the entropy-generation rate when sources are 
translated. 


A. A body in an external field—extension of the d’Alembert paradox 


We saw that when V<p vanishes at infinity the total force on any bound system vanishes. When we have V(p —> —go 
at infinity, describing a constant external field in which the system is immersed, it can be shown that even in the 
nonlinear case the force on the whole syste m is Qgo, where Q is t he total charge. This can be seen by using the 
asymptotic form of the field given below eq.( 2.24 ) in expression (^) for F on a surface going to infinity. In the linear 
case this result is trivial since the total force is the sum of the total mutual forces of the charges, which vanishes, and 
the sum of the external forces f pgg d’^r = Qgo- It follows from this, for example, that the force on an arbitrary 
/i inclusion in a constant external field vanishes since, as we saw above, its total effective charge vanishes. The 
application of this result to the compressible flow problem constitute a generalization of the well-known d’Alembert 
paradox in fluid mechanics to the effect that a single, static obstacle in a non-viscous, incompressible, potential flow 
that is uniform at infinity is subject to no force. We now see, more generally, that the paradox applies as well to the 
nonlinear (compressible but subsonic) case, and to any single body made of fluid sources (or sinks), rigid obstacles, 
regions where body forces apply, or any other configuration that can be replaced by an effective charge distribution 
with vanishing total charge. 


B. Point bodies 


When a body is very small compared with the typical scale of the system we can approximate it by a point charge 
(PC) q at position r^, with density p{r) = qS^{r — Tq). We may view the body as the limit of some finite-size charge 


p(r) = g lim A 
A^O 




(3.5) 


where p is some smooth, finite charge distribution normalized to / /5(R)d^i? = 1. Not all problems admit point 
charges. Applying Gauss’s theorem to eq. (|l.l| ) we see that wp,{w) diverges as near a point charge, so such 

charges are not admitted in theories with wp{w) bound from above. This is the case for the volume-minimization 
problem, where wp,(w) < 1. It is also the case in the flow problem with > 0, where subsonicity, and hence ellipticity, 
is lost near a point source. For “flow” problems with < 0, ellipticity is maintained at all values of Vtp, and point 
sources are not objectionable. 

The concept of a PC is useful only if the field everywhere in a system containing a PC is independent of the 
particular choice of the structure function p of the PC, so that it enters only though its total charge. I was not able 
to prove that this always holds. One possible route to proving this is to show that an infinitesimal change <5/5 in 
the structure function, which does not change the total charge, f Sp(r) d^r = 0, prod uces a n everywhere-vanishing 
increment in the potential. The potential increment C is a solution of the linear equation (2.13) (where the background 
field depends on A) with 6p = X~^6p{r/X), in the limit A —> 0. In this limit all the moments of Sp vanish, which might 
tell us that C vanishes: Start with the charge distribution p\x' = X~^p{r/X) + X’~^5p{v/X') in place of the point PC. 
Let ip\ be the solution of the problem for 6p = 0, and -f Caa' the solution for pw. In the limit A —> 0, goes 
to the solution for the PC having /3 as a structure function. In the limit A = A' ^ 0, </3 a + Caa' goes to the required 
solution for a PC with structure function p + 5p. We want to show that in this last limit Caa' ^ 0. Instead of going 
to the limit A = X' = 0 along the A = A' line we first take the limit A' ^ 0 for finite A and then take the limit A —> 0. 
Caa' solves eq.( 2.13 ) with the left-hand side depending on A, and the right-hand side source being X'~^5p{v/X'). As 
all the moments of this source vanish in the limit A^ ^ 0, I conclude that ^aa' —> 0 in the limit. Taking now the limit 
A —> 0 we are left with (j = 0 in the limit. The remaining loophole concerns the equality of the two limits. 
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Be this as it may, the results concerning PCs are only valid when t he h eld does not depend on the structure 
function. Employing the expression for the force as a surface integral eq.(^) we see that the force on a PC is then 
also independent of the choice of p. 

Zero-size bodies with higher multipoles can, of course, also be constructed by considering a multiple limit with 
point charges of infinite charges at zero distances. Such bodies with pure multipole charge distribution do not have 
the special role they have in the linear case. They do not, in general, produce unique-multipole fields, and, unlike the 
pure-charge case the force on them depends on details of their structure, not only on the components of the multipole. 


C. Test bodies 

Consider a sub-system of charge distribution p^- In the limit Pb ^ 0 it can be considered as a test body: its 
contribution to the full potential can be neglected in expression ( |3.3| ) for the force on itself, and we have F = — f p^Vp, 
where p is the potential determined by the rest of the system. A point charge has an inhnite density. But, it may 
still be considered a test charge (the force on which is F = —qV(p) provided q is small enough that a surface can be 
drawn around it such that (1) the surface is small compared with the scale over which the held varies appreciably, (2) 
it is far enough from the charge that the latter’s contribution on the surface is small compared with the held due to 
the rest of the system alone (which is then approximately constant on the surface). 


D. Attraction or repulsion? 


I now show that it is still correct in the nonlinear case that two equal point charges attract, and two opposite 
charges repel each other (for G > 0; and vice versa for G < 0). Take, more generally, two hnite, disjoint bodies with 
charge distributions that are the mirror images of each other about some {D — l)-dimensional hyperplane. The force 
on each can be calculated by taking the symmetry hyperplane as the integration surface in eq.(3.4) (the integration 
on the hemisphere at infinity vanis hes with our choice .F(O) = 0). Since p is symmetric about the hyperplane, Vp is 
in the plane. We thus get from eq.(3.4) 


^ Qq 
2aoG 


T ds. 


(3.6) 


which is attractive for G > 0. If one body is the negative-charge reflection of the other, p is antisymmetric, and Vtp 


is perpendicular to the hyperplane so that "P • ds oc F{1 — 2iF)ds, and from eq.(3.4) 


F = 

J- -7) 


2a„G 


T{1 - 2T)ds, 


(3.7) 


which is repulsive for G > 0 since JF > 1/2 when Vp ^ 0. 

Consider, more generally, two disjoint bodies Bi defined by charge distributions pi > 0 in the non-overlapping 
volumes vi, i=l,2. It is meaningful to ask whether they attract or repel each other only if they can be separated by 
some {D — l)-dimensional hyperplane. We have attraction if the force on Bi crosses any such separating hyperplane 
from the side of Bi to that of B 2 (see Fig. 1). I conjecture without a general proof that indeed this is always the 
case, and that if the two bodies are oppositely charged, i.e. if, say, p 2 < 0, they always repel each other. (The sign of 
the charge within each body must be uniform.) It follows from the result proved in Appendix C that the conjecture 
holds when one of the bodies, say i?i, is spherically symmetric with a density profile decreasing from the center out 
(I thank Shoshana Kamin for discussions leading to this proof). In fact, in this case the statement is stronger: for 
any hyperplane H* through the center of Bi with B2 wholly to its one side (side 2) the force on Bi points from side 
1 to side 2 (see Fig. 2). 

We learn from the above result that the force on a spherical body with a decreasing density profile is always within 
the convex closure of the cone defined by its center and the other body (the envelope of all the planes through the 
point that are tangent to the body). 

Because a point charge may be considered as the limit of such a spherical body, we deduce that the force on a point 
charge in the presence of an extended body of uniform-sign charge is always within the cone from the point to the 
convex closure of the body. In particular, two point charges of the same sign always attract each other. 

All the above generalizes, mutatis mutandis, to oppositely charged bodies. 
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FIG. 1. The force F on body 1 crosses every separating plane from the side of 1 to side of 2 



FIG. 2. The setup for demonstrating attraction of like charges when one of the bodies is spherical with a decreasing density 
profile 



FIG. 3. A conhguration depicting the push-pull conjecture: for any pair of parallel separating planes the force F on Bi 
crosses the planes into side 2 
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The two body case may be generalized to a push-pull conjecture concerning three bodies (as in Fig. 3) with 
Pi, P 2 > 0, p 3 < 0: for every two parallel hyperplanes separating the three bodies the force F acting on Bi crosses 
the planes into the side of Bi. This seems to be a good way to summarize the concept of attraction-repulsion of 
bodies with constant-sign charges as all else regarding the question follows from it. 

I have not been able to prove this push-pull conjecture in all generality. It is proved in Appendix C when body 1 
is spherically symmetric with radially decreasing density; so, in particular, the conjecture holds when Bi is a point 
charge. The conjecture is also proved in the one-dimensional case (see below). It also holds, generally, when body 1 
is a collection of test charges. It then follows from the comparison principle discussed in Appendix B. In this case the 
force on body 1 is F = — pV(pi 2 , where ipi 2 is the potential produced by bodies 2 and 3 alone. The comparison 
principle tell us that on any plane separating these bodies —V(pi 2 crosses from the side of 3 to that of 2 so Fi does so 
as well. Similarly, the conjecture holds when bodies 2 and 3 are made of test particles. The force on body 1 can then 
be written as F = pV'Pi, where the integration is performed over bodies 1 and 2, and (pi is produced by body 1 
alone. By the corollary of the comparison principle, W(pi crosses from the side of 3 to that of 2 in the volume of body 
2 , and the opposite in the volume of body 3. Thus, the integrand, and hence F, crosses from 3 to 2. 

Take a system containing charges of a uniform sign, say p > 0, with § a convex surface surrounding all the charges. 
It follows that on S', Vtp points outward (inward when p < 0). So, for example, in theories where /i < 1 we have 
for the total charge, Q, the inequality Q = {aoG)~^ pVp • ds < (ajyG)~^ Vp • ds = Q*. Q* is the charge that 
would be deduced from the held Vp in a linear theory. 


E. Forces in the one-dimensional case 


The one-dimensional case can be solved in closed form once the boundary conditions are hxed. If we require that 
V<p(oo) = —V(p(—oo) so as to nullify the force on an isolated system, then Vp ^ const, at oo. The force, F, on 
an arbitrary charged body of total charge q, in the presence of a charge distribution that does not overlap with it 
depends only on q and on the difference, Q, between the total charges to the right and to the left of the body. 


F = Fiq,Q) = 


^[x(w+)-x(w )], 


where x(w) is dehned in eq.(^), and 


(3.8) 


=s(G)z/ ^ (3.9) 

We see that F is invariant under translations of the body, as long as it does not cross other charges, and that 

F{-q,Q) = -F{q,Q)=F{q,-Q), F{Q,q) = F{q,Q). (3.10) 

These are peculiarities of the one-dimensional case. In higher dimensions the magnitude of the force on a charge does 
not, in general, remain invariant when the sign of the charge is reversed (all others kept intact). 

Since x' > 0 for w > 0, x is increasing everywhere. Thus, F does not vanish unless w'^ = w~, i.e., unless g = 0, or 
Q = 0. The push-pull conjecture hold here as is easily seen. 


IV. FORCES ON POINT CHARGES 
A. The limit of a very large point charge 

Consider a charge distribution p and a point charge q that does not overlap with it. Assume also that q is much larger 
than any partial charge that makes p up (i.e. p d^r q for any volume v). We can then consider p to be a collection 
of test charges relative to q. From momentum conservation, the force Flj(R) acting on q at position R is opposite the 
force Fp acting on the distribution p. This latter is given in the test-particles limit hy Fp = — f d^rp(r)V(pq, where 
ifq is the held produced by q alone, which is straightforwardly gotten through Gauss’s theorem to give 

F,(R) = ,(,G)a„ I d-ri,{v)u-' (4.1) 
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with v{w) = w^{w). This force is derivable from an effective potential Eqi 


Eq{R)= I d^rp{r)gq{\r-R\). (4.2) 

where the effective Green’s function, Qg, satisfies 

^rGgir) = s{qG)aQV~^{z)n, (4.3) 


with z = \qG\/ao\r\’^ \ and n = |^. 

The effective potential is then a linear functional of the density p. It satisfies the Poisson equation 

^Eq{R) = j d^rp{v)/\gq{\v-R\), (4.4) 

with 

A5,(r) = s{qG){D - l)ao|rri^^^^, (4.5) 


where w = 

The above results apply also to any finite spherical body of very large total charge, and are straightforwardly 
extended to the case of an arbitrary body whose ip field is known. 


B. The two-body force 

For two point charges ( 7 i < 92 , a distance ^ apart, the vanishing of the total force and moment tell us that the 
forces on the two charges are opposite each other, and lie along the connecting line; write its magnitude /( 9 i, 92 ,f')- 
We have seen above that / has the sign of G 9192 (/ is positive for attraction). Because £ is the only length scale, we 
can reduce the number of variables to two independent, dimensionless variables constructed from 91 , 92 , £, and oq; for 
example, —1 < 77 = 91/92 < 1 and z = |G 92 t/ao£°“b Since ao 9 i has dimensions of force, we can write, for example, 

/ = s{Gq 2 )aoqJ{r], z). (4.6) 

Our earlier discussion implies some constraints of /. For instance, when one charge is much smaller than the other 
ryl ^ 1 , the test-charge result tells us that 

/(|77 |«l,z) = z7-i(z), (4.7) 

to lowest order in 77 [z = v{w) = wp{w)]. 

For a power-law medium with p{w) = one deduces from scaling properties of the field equation that the 
^-dependence of the force is where 7 = 1/(1 -I- /3). So here 

/ = (4.8) 

From eq. © C( 0 ) = 1 . Since, in general, C(~l) 7 ^ C(l)) the forces for equal and for opposite charges are not of the 
same magnitude, in contradistinction with the linear case. For the special case (3 = D — 2, the two-body force was 
found in closed form 

/(9i,9.,^) = s{G)jd-^\Ga^S~\k^ + - kil" - k.l") (4.9) 

[d = D/{D — 1) = 1 -I- 7]. (The forces in a three-point-charge system of zero total charge were also derived in [^.) 
For two equal charges 9 i = 92 = 9 , / = 2s(G)£“^fi“^|G|‘’*“^| 9 |'’*( 2 ‘^“^ — 1), while for opposite charges: 91 = —92 = 9 , 
/ = — 2 s(G)£“^(i“^|G|‘^“^| 9 |'’*. The two are equal in magnitude only in the (linear) two-dimensional case. Interestingly, 
in the li mit o f large dimension, where d —>■ 1, the two-body force for two charges of the same sign vanishes as D~^ 
[from eq. (|4.9[ )], while the force for opposite charges does not. This can be generalized: for a given configuration of 
N point charges of the same sign the force on each becomes smaller as D~^ in the limit of large D (letting also the 
theory’s power increase, (3 = D — 2). 
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V. SOME EXAMPLES OF APPLICATIONS 


I next discuss a number of potential applications, some of which I alluded to earlier. I concentrate on two of the 
physical problems listed in the introduction, viz stationary (potential) flows of barotropic compressible fluids, and 
media with field dependent dielectric, or diamagnetic, constants. 

Concerning the flow problem, we saw that the d’Alembert paradox can be extended to the compressible case. So, 
the force acting on a body of total effective charge in a flow of constant asymptotic speed vanishes. This implies the 
vanishing of the force on a body of rigid walls standing in such flow, on a system made of sinks and sources with 
vanishing total out-flux. If the fluid is electrically charged-weakly, so that the fluid does not self-interact-then the 
flow is modified by the presence of a region with an electric field. The effect of such a region can be described by 
an effective sink-source distribution whose total out-flux vanishes. So, the net force on such a region due to the fluid 
motion vanishes. We also learned from the discussion of the sign of the forces that in such flows two sinks, or two 
sources always attract each other while a source and a sink repel each other. 

For the nonlinear dielectric, the generalized d’Alembert paradox says that, as in the linear case, the force on an 
arbitrary charge distribution, e^, of vanishing net charge in a constant external field, E, vanishes, even though the 
force on an individual component of charge is not eiE. From this we learn that, for example, the force on a dielectric 
inclusion, or an equipotential body of zero net charge, vanishes in a constant external field. We also deduce that the 
force vanishes on a superconducting inclusion in a nonlinear magnetic medium in a constant magnetic field. 
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A. INTEGRAL RELATIONS 


We employ the techn ique described in to derive certain useful virial relations directly from the action by 
substituting in eq.(2.I5) for 5S various choices of 5ip. If either /3o < D — 2, or /?(, = D — 2 and (5 = 0, the potential 
vanishes at infinity. We can then obtain one relation by taking 5{p = with e infinitesimal. The vanishing of (f at 
infinity leads to the vanishing of 5S. But 5S can also be calculated directly to yield a virial relation 


p(p 


1 


O-dG , 


{y(pYp.{\V(p\/ao) d°r = 0 . 


(A.l) 


(This can also be derived by multiplying p by the expression for p from the field equation, and integrating by parts.) 

In the relation above (and below) p should be understood to include true sources as well as all the effective sources 
replacing boundary conditions on closed surfaces in the system. 

Another useful relation, satisfied by solutions of the field equation, may be obtained by considering variations 
of S produced by dilations of space coordinates: p{v) —> (/?(Ar). For an infinitesimal dilation A = 1 -I- e, we have 
6p = er ■ Vp. Now, dS can also be calculated directly by substituting p{y) — > :p(Ar) in S'; then taking the derivative 
with respect to A at A = I to obtain for the finite volume 


lx 


= — pr ■ Vp d°r + 


A=1 


2anG 


[DT-2 


(ypY 


p] d^r - 


2anG 


Tr ■ ds. 


(A.2) 


So, comparing with expression (p . 1 5D for the variation, we get an expression for the virial V in the volume V. 


V = / pr ■ Vp d^r = 


2a„G . 


T{D -2T)d^r+ rP- ds. 


(A.3) 


This can also be obtained by substituting in the definition of V the expression for p from the field equation, and 
integration by parts. The value of V is independent of th e cho ice of origin due to momentum conservation. 

Taking the surface to infinity, the surface term in eq.(A.3) can be evaluated. At this point I restrict myself to 
D >2; the one-dimensional case is exactly solvable and is treated below. For Pq < D — 2 the surface integral vanishes 
at infinity, and 
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(A.4) 


For Po = D — 2 the surface integral at infinity converges to yield 

V = [ T{D - 2T) d°r + dT^G-pCQ^^, (A.5) 

where d = D/{D — 1), and G = Goq®. The conformally invariant theories discussed in have T = D/2 and so 
V = d~^G~^\GQ\'^. This has been put to extensive use in |^. 

Another integral of interest is 

a = f p[r(r ■ Vip) — d°r. (A.6) 

Jv ^ 

By substituting p from the field equation and integrating by parts we arrive at the expression 
-> rP r - /■ 

U = -—^ / T{D — 2T)v d'^r — -—^ — 2n ® n) + ^[4(n • e)n 0 e — 2e 0 e]} • ds (A.7) 

2q!dG Jy 4o;j:,G Jj, 

(n = r/|r|, e = Vv?/|Vv?|). 

Note that U depends, in general, on the choice of origin. Using the vanishing of the total force and total moment 
we see that if the origin is shifted by —a, U is changed by Va. 


B. THE COMPARISON PRINCIPLE AND SOME CONSEQUENCES 


First I show that if (pi, ip 2 are continuous functions that solve our equation for densities pi > p2 in a volume V of 
boundary S, and ipi < ip 2 on S, then ipi < ip 2 everywhere in V. This is known as a comparison principle for solution 
of elliptic equations (e.g. |^). I give here a proof that applies specifically to the form of our equation and is thus more 
elementary than the proofs found in the literature. Start with the identity 


{<Pi - ‘P 2 ){pi - P 2 ) d°r cx {<pi - (^2)V • {piV(pi - P 2 ^<P 2 ) d°r = 


(B.l) 



tp2){Pl^‘Pl - P2^‘P2) ■ ds 



V(P2) • - P2^P2) d°r 


(B.2) 


[where pi = /r(|V(/?i|)], gotten by using the expression for pi from the field equation, use of Gauss’s theorem, and 
integration by parts. I want to show that the region v in U, in which pi > p 2 is empty. If there is even one point 
where pi > p 2 , then, by continuity of pi — p 2 , v must contain a whole (non-zero-measure) neighborhood. Apply 
identity |B.2| to this volume v. On its boundary a we have pi = p 2 (from continuity of i^^), whether a has overlap with 
E, or is completely interior to V. Thus the first term on the right-hand-side of p3.2| vanishes. It can be shown (see |^) 
that, in light of the ellipticity condition, the integrand in the second term is non-negative, and vanishes only where 
Vpi = Vp 2 - This is because when wp{w) is non-decreasing (a — b) • [/x(|a|)a — ^(|b|)b] > 0 for any two vectors a, b, 
and vanishes only for a = b. However, by the assumptions, the left-hand side is non-negative and so the right-hand 
side must vanish, and hence V(^i = V(^2 in v. Thus, pi — p 2 is constant in the whole region where it is positive. 
This, however, contradicts the assumption that pi < p 2 on If with the continuity of the potentials. 

I now apply the theorem to the following useful configuration. Let p{xi, ...,Xo) be the solution for a source 
distribution that satisfies p{xi, ...,Xo) > pi—xi, ...,Xd)- The boundary condition at infinity is p(v) sir). Look on 
the two halves of the solution in the two half-spaces separated by the xi = 0 hyperplane {H*) as two solutions of 
the field equation in a half-space Xi > 0; so, pi(xi, ..., Xd) = p(xi, ..., Xd), and p 2 {xi, ..., Xo) = piy—xi, ..., Xo)- We 
have (^2(0, a;2,..., a;^,) = (/5i(0, Xa,..., a;^). So, pi and p 2 are solutions with the same boundary values (they also have 
the same boundary condition at 00). The comparison principle then tells us that pi(xi, ...^Xd) < p 2 {xi, ...,Xd), or 
p{xi, .■.,Xd) < ‘fi—xi, ...,Xd)- In particular, J^(a;i = 0) < 0. H* can be any plane separating a source distribution 
pi > 0 from p 2 < 0. 
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Let /9 > 0 in a volume v be the whole the source distribution, and A its convex closure {A is the smallest convex 
volume containing all the points where p ^ 0). Then the gradient of the potential at any point r outside A points 
away from A (because it points away from the side of A on any plane separating r from A). 

We can deduce from this that for a system of two point charges of the same sign the tangent to the field lines alway 
crosses the line connecting the charges between the two. For oppositely charges points the tangent crosses outside 
this line segment. This generalizes the usual situation familiar from the linear case where it follows simply from the 
vector addition of the two forces due to the point charges. 


C. PROOF OF THE PUSH-PULL CONJECTURE FOR A SPHERICAL BODY 


Consider a system of three disjoint bodies: Bi, which is spherically symmetric with pi > 0 decreasing from the 
center out; B 2 with p 2 > 0 ; and S 3 with pa < 0 . H* is a plane through the center of Bi that separates B 2 from S 3 . 
I show that the force F on Si crosses H* from side 3 to side 2. The push-pull conjecture for this special case is a 
weaker statement and follows as a corollary. 

Choose the coordinates such that H* is the xi = 0 plane, with S 2 on the xi < 0 side. We then have from the 


previous appendix that ip{—Xi^ ■■■,Xo) < ...,Xo) for xi > 0. Using the second expression in eq.(^) to calculate 

the xi component of the force on Si, and employing the latter’s symmetry, we have 


Fl= d°r [ip{xi, ...,Xd) - pi-Xi, ...,XD)]da:iP <0. 

J a:;i>0 

where use was made of the fact that dx^p < 0 for xi > 0 , as p is monotonically decreasing with radius. 
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